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Abstract. We consider a magnetic Schrodinger operator H'^ , depending on 
the semiclassical parameter /i > 0, on a two-dimensional Riemannian mani- 
fold. We assume that there is no electric field. We suppose that the minimal 
value feo of the magnetic field b is strictly positive, and there exists a unique 
minimum point of b, which is non-degenerate. The main result of the paper is 
a complete asymptotic expansion for the low-lying eigenvalues of the operator 
in the semiclassical limit. We also apply these results to prove the exis- 
tence of an arbitrary large number of spectral gaps in the semiclassical limit 
in the corresponding periodic setting. 



1. Preliminaries and main results 

Let M be a compact oriented manifold of dimension n>2 (possibly with bound- 
ary). Let g be a Riemannian metric and B a real- valued closed 2-form on M. As- 
sume that B is exact and choose a real- valued 1-form A on M such that dA = B. 
Thus, one has a natural mapping 

u i-^ ih du + Au 

from C^{M) to the space fll{M) of smooth, compactly supported one-forms on 
M. The Riemannian metric allows to define scalar products in these spaces and 
consider the adjoint operator 

{ihd + A)* : nl{M) ^ C^{M). 

A Schrodinger operator with magnetic potential A is defined by the formula 

ij'' = {:ihd + Ay(ihd + A). 

Here /i > is a semiclassical parameter. If M has non-empty boundary, we will 
assume that the operator H'^ satisfies the Dirichlet boundary conditions. 

We are interested in semiclassical asymptotics of the low-lying eigenvalues of the 
operator . This problem was studied in[3l[8l[l0l[ITl[l2l[l3l[l4l[23l[25l[26l[2Zl 
nil [201 [SI] (see [l[n] for surveys). 



2000 Mathematics Subject Classification. 35P20, 35J10, 47A75, 58J50, 81Q10. 
Key words and phrases. Spectral theory, Schrodinger operators, magnetic fields, eigenvalue 
asymptotics. 

Y.K. is partially supported by the Russian Foundation of Basic Research (grant 09-01-00389). 

1 



2 



BERNARD HELFFER AND YURI A. KORDYUKOV 



In this paper, we study the problem in a particular situation. We come back to 
the case considered in [11]. We suppose that M is two-dimensional. Then we can 
write B = bdxg, where h € C°°{M) and dxg is the Riemannian volume form. Let 

60 = min bix). 

xeM 

We furthermore assume that: 

(1) &o > 0; 

(2) there exist a unique point xo, which belongs to the interior of M, fc G N 
and C > such that for all x in some neighborhood of xq the estimates 
hold: 

C-^ d{x, xof < b{x) -bo<C d{x, xof ■ 

Denote 

1 V/' . /I 



a — Tt ^-Hess 6(a::o)J , d — det |^-Hess 6(a;o)J . 

Denote by Ao(i/'') < Xi{H^) < A2(iJ'') < ... the eigenvalues of the operator 
H'' in L^{M). 

Theorem 1.1. Under current assumptions, for any natural j, there exist Cj > 
and hj > such that, for any h € (0, hj], 



hbo + 



bo ^ ^ 2bo 



Cjh^^/^ < Aj(i7^) < hbo + h^ 



2d^. 

bo ^^2bo 



In particular, we have lower and upper bounds for the groundstate energy 

XoiH'^y. 

2 2 

hbo + h^-^-Coh^^/'' <Xo{H'') <hbo + h^^+Coh^^^, /iG (0,/io]. 
ZbQ Zbo 

and the asymptotics of the splitting between the groundstate energy and the first 
excited state : 

oo 

The previous statement can be completed in the following way. 

Theorem 1.2. Under current assumptions, for any natural j , there exists a se- 
quence {aj^e)eeN, and for any N , there exist Cj.N > and hj,N > such that, for 
any h G (0, ft-j,Ar], 

N 

(1.1) \Xj{H'') ~ hJ2o^j,ih^\ < Cj,Nh^ , 

e=o 

with aj^o = bo, ajs = 0, 0^,2 = ^f^-J + ^ • 

This theorem improves the result of [llj which only gives a two-terms asymp- 
totics for the ground state energy in the flat case. 

The scheme of the proof is to first prove the weak version, given by Theorem 1 1.1[ 
with N = 2, permitting to determine jo disjoint intervals in which the first jo 
eigenvalues are localized, for h small enough, and then to determine a complete 
expansion of each eigenvalue lying in a given interval. 

The paper is organized as follows. In Section [2] we construct approximate eigen- 
functions of the operator with any order of precision. This allows us to prove 
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accurate upper bounds for the jth eigenvalue of H^. In Section [3] we prove a lower 
bound for the jth eigenvalue of and complete the proofs of Theorems l 1 . ll and ll .21 
In Section m we consider the case when the magnetic field is periodic. We combine 
the construction of approximate eigenfunctions given in Section [2] with the results 
of [6] to prove the existence of arbitrary large number of gaps in the spectrum of 
the periodic operator H'^ in the semiclassical limit. 

We wish to thank the Erwin Schrodinger Institute in Vienna and the organizers 
of the conference in honor of 65th birthday of Mikhail Shubin "Spectral Theory 
and Geometric Analysis" in Boston for their hospitality and support. 

2. Upper bounds 

2.1. Approximate eigenfunctions: main result. The purpose of this section 
is to prove the following accurate upper bound for the eigenvalues of the operator 

Theorem 2.1. Under current assumptions, for any natural j and k, there exists 
a sequence {fJ.j.k.e)iefi with 

ti-j.k.o = (2fc + l)&o, t^j.k.i = 0, 

and 

f^j,k,2 = {2j + l)(2fc + 1)^ + (2fc2 + 2k + 1)-^ + i(fc' + k)R{xo) , 

Oo /Do Z 

where R is the scalar curvature, and 

t = TrQHess5(xo)^ , 
and for any N, there exist <j)'^i.j^i € C°°{M), Cjk,N > and hjk^N > such that 

(2-1) (0jifei7V>i2fe2Ar) = Sjl32Sklk2 + Cji j2,fei,fe2(M 

and, for any h G (0, hjk^N], 

(2.2) \\H''<f>%N - f^%N<t>%N\\ < QkNh'^U^^l 

where 

N 

1^%N = hJ2lJ-j,k,eh^- 

Since the operator H'^ is self-adjoint, using Spectral Theorem, we immediately 
deduce the existence of eigenvalues near the points /i^^, . 

Corollary 2.2. For any natural j , k and N , there exist Cjk,N > and hjk^N > 
such that, for any h € (0, hjk,N), 

dist(Ai^\^,Spec(i7'*)) < Cjk,Nh^. 

Remark 2.3. The low-lying eigenvalues of the operator as /i — ?> 0, are obtained 
by taking fc = in Theorem 12.11 Therefore, as an immediate consequence of 
Theorem 1 2. 1[ we deduce that, for any natural j and N , there exists hj^M > such 
that, for any h G (0,/ij,Ar], we have 

— t^%N + CjQ,Nh 2 . 
In particular, this implies the upper bound in Theorem ll.il 
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Remark 2.4. Our interest in the case of arbitrary k in Theorem 12.11 is motivated 
by its importance for proving the existence of gaps in the spectrum of the operator 
in the semiclassical hmit. This wih be discussed in Section H) 



Proof of Theorem \2.1\ The proof is long, so we will split it in different steps in the 
next subsections. □ 

2.2. Expanding operators in fractional powers of h. The approximate eigen- 
functions G C°°(M), which we are going to construct, will be supported in 
a small neighborhood of xq. So, in a neighborhood of xq, we will consider some 
special local coordinate system with coordinates (x, y) such that xq corresponds 
to (0,0). We will only apply our operator on functions which are a product of 
cut-off functions with functions of the form of linear combinations of terms like 
h'^w{h~^^^x, h^^^^y) with w in 5(]R^). These functions are consequently 0{h°°) 
outside a fixed neighborhood of (0,0). We will start by doing the computations for- 
mally in the sense that everything is determined modulo 0{h°°), and any smooth 
function will be replaced by its Taylor's expansion. It is then easy to construct non 
formal approximate eigenfunctions. 

First, we recall that in local coordinates X — [X^ ^X"^) — (x, y) on M the 1-form 
A is written as 

A = Ai{X) dXi + A2{X) dX2 , 
the matrix of the Riemannian metric g as 

9(X) = {gjdX))i<j.e<2 

and its inverse as 

g{X)-' = {g^\X)),<,,,<2. 
Denote 1 5 (-'4^)1 = det{g{X)). Then the magnetic field B is given by 

B = BdxAdy, B^^-^, 
ox ay 

and 

B = b^\. 
Finally, the operator has the form 



where 

or, equivalently. 



<Q,^<2 



ih^+A^{X), a = 1,2, 



l<a,P<2 l<a<2 

where, for a = 1, 2, 



We will consider normal local coordinates near xq such that xq corresponds to 
(0, 0) and, in a neighborhood of xq, 

h{X) ^bo + aix^ + + Oi\X\^). 
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Thus, we have 

a = (ai)i/2 + (/3^)i/2^ d = t = ai + Pi. 

By well-known properties of normal coordinates we have 

gn{X) = l + 0{\Xf), gi2{X) = 0{\Xf), g22{X) = 1 + 0{\Xf). 
Moreover (see, for instance, [1] Proposition 1.28]), we have 

kl 

where Riju is the Riemann curvature tensor. Therefore, Taylor's expansion of g"^ 
has the form 



oo 



(2.4) g'^^iX) = S'^P + J29'^^^iX), 

k=2 

where 

9li^{x) = \Y.^-m{^^)x''xK 

kl 

In the two-dimensional case, due to its symmetries, the Riemann curvature tensor 
is determined by a single component 

^1212 = — -R2112 — -R2121 = — -Rl221- 
The other components equal zero. We have 

2i?i2i2 = ^(511522 - 512): 
where R is the scalar curvature. So we have 

i?i2i2(a;o) = ^R{xo)- 



Thus we have 

(2.5) gll){X) = jR{xo)v\ gl^){X) = -^R{xo)xy, gf^^iX) = ^R{xo)x' 



We also have 



1 . . 9 1 



(2.6) VWni = 1 - -Rixo)x' - j^R{xo)y^ + 0(|X|3). 

Let us write Taylor's expansion of in the form 

00 

(2.7) r"(x) = ^rf,)(x), a = 1,2. 

fc=0 

Using ()2.5p and ()2.6p . one can show that 

r^o)(x) = o, a = 1,2, 

and 

(2.8) Tl^iX) = -^R{xo)x , rl^iX) ^ -^R{xo)y . 
If we write — B(x,y) dx dy then 



B{X) = b{X)MX)\ - 60 + ax^ + py^ + 0(|X|3), 

where 

1 1 

ai^a + —boRixo) > 0, /3i = /3 + — 6o^(2;o) > 0. 



6 



BERNARD HELFFER AND YURI A. KORDYUKOV 



We can also choose a gauge A such that 

Ai{X) = and A2{X) = 60a; + ^x^ + I3xi/ + 0{\X\'^). 
We expand A2 into the Taylor series: 

00 

A2{X) ^ box + J2Sj{x,y), 

3=3 

with 

In particular, we have 

S3{x,y) = ^x'^+pxy^ 

Next we move the operator H'^ into the Hilbert space L^(K") equipped with the 
Euclidean inner product, considering the operator 

H'^ = \giX)\'^'H,MX)\-'/^ 

= ^ g"P{x)w'y'^ + ih ^ r-v^, 

l<Q,/3<2 1<Q<2 

where, for a = 1,2, 

(2.9) = |ff(X)|VX|.9(X)rV4 = + ^^hR{xo)X- + 0{h\X\')). 

Now we use the scaling x = h^^^xi, y = h^^^yi and expand the resulting operator 
H'^{xi,yi, Dxi, Dy^) into a formal power series of h^^"^. By (|2.9p . we have 

= h^'\^D^, + ^ihR{xo)xi + 0{h^^^)), 

= + boxi + -^thR{xo)yi + hSsix^^yi) + 0{h"^)). 

From ([2?4)) and ([2J)) . we get 

oo 

r"(xi,yi) =/^i/2pa^(a;i,yi) + ^/i'=/2rffc)(xi,yi), a = 1,2. 

Using these expansions, one can check that the operator has the form 
H''{xi,yi,D^„Dy,) = hQ''{xi,yi,D.^„DyJ, 

with 

CXD 

Q^{xi,yi,D,„Dy,)^Yl h'^'^Qk ixi,yi, Dy,), 

where 

QQ{xi,yi,D^^,Dy^) = Dl^ + [Dy^ - b^xif, 
Qi(xi,yi,D^^,Dy^) = 0, 
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Q2{xi,yi,Dx^,DyJ = - —iR{xo)ixiDa:^ + D^^xi) 



and 

1 

- ■Y^'i'R{xo){yi{Dy^ - boxi) + {Dy^ - boxi)yi) 

- {{Dy^ - boXi)S3{xi,yi) + Ss{xi,yi){Dy^ - boxi)) 

+ 9l2)i^i,yi) D^, {Dy^ - boxi) 
+ 9{2)ixi,yi) {Dy, - 6oa;i)Dxi 
+ 9'(2)i^i,yi){Dy, -boXif 

- irf^^{xi,yi)(Dy^ -boxi). 

Using the Fourier transform in yi, we can write the operator Q'^ as 

00 

fe=0 

where 

Qo{xu-D^, D,„0 = Dl^ + - box^f, 
Qi{xi,-D^,D^„O=0, 

and 

Q2{xi,-D^,D^^,^) =- ^iR{xo){2xiD^^ -i) 

+ ^tR{xo){2{^-baXi)D^-i) 

- ((C - box,)S3{x,, -D^) + S3ixi,^D^)i^ - box,)) 
+ g\l^{x,,^D^)Dl^ 

+ 5(2) (^1' ~^«) -^^1 ~ ^0^1) 

+ gl^){xu-D^){^-box,f 
-iT\^){xi,-D^)D^, 
-iTf^^{xu-D^){^-boXi). 
A further translation X2 = Xi — ^ gives 

Hf^ = hT^{x2,^,D,„D^), 

where 

T''{X2,^,D,„D^) = (^X2 + ^,-D^ + 1d,„£),„^^ . 
We have (denoting w = {x2,0) 

00 

T^{w, D^) = To{x2, D^,) + hT2{w, D^) + ^ h^'%{w, D^) 

j=3 



8 



BERNARD HELFFER AND YURI A. KORDYUKOV 



with 
and 

T2{w, DJ = - ^iR{xo){2 (^X2 + D,, - i) 

+ ^iR{xo){2boX2 (^-D^ + l^D,2^ + i) 

+ {boX2Ss{x2,^,D^^,D^) + Ss{x2,^,D^^,D^)boX2) 
-.11 r„ , ^ n , 1 n \n2 



X2 



+ 9\l^{x2 + ^^,-D^ + yD^.)D 

- 5(2)(^2 + ^, -D^ + ^-Dxa) {boX2) 

- 9{i){x2 + -D^ + ^D^^) boX2 -Dx2 

+ 5(2) (^2 + -D^ + ^D,,)blxl 

Oo Do 

- iTl^ix2 + + ^^-2)^x2 
+ iTl)ix2 + f,-D^ + ^D^,)boX2, 



where 



Ss{x2,^, D^,,D{) = Si[x2 + Y^,-D^ + 



The operator ^3 has the following form: 

Si{x2,^, D,„D^) = X2L{i, D^) + Mo(e, D^) + M^{x2. D,,) 



+ M2{x2,D^„D^) + M3 {X2 ,D^„D^) + M4{X2,^,D^,), 



where 



a 



"0 

Mo(^,^,) = |3^^ + £^i^|, 
Mr{x2,D^^)='^xl + ^X2Dl^, 



M2{X2,D,,,D^) = -2^D,,CD^, 
"0 

M3ix2,D^„D^) = -2^X2D^,D^, 
oq 

M4x2,^,D,,) = ^xU + l^D^ 



bo bl 



X2 ■ 



So we get 



T2{w,Dyj) = - ^iR{xo)-^Dx^ - ^iR{xo)boX2D^ 

+ 2boxlL{^, D^) + bo {x2M{w, D^) + M{w, D^)x2)) 
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where 



We have 



Sp(To(a;2, D^,)) = {fik = {2k + l)6o : e N}. 
The eigenfunction of To(a;2, -DX2) associated to the eigenvalue fih is 

M^2)=7:-'/X/'H,{bl/'x2)e-''°-i/^ 
where Hk is the Hermite polynomial: 

The norm of V'fe in L^{^,dx) equals the norm of H/. in L^{M.,e~^ dx), which is 
given by 



\\Hk\\ = \l2kk\Vk. 



2.3. Construction of approximate eigenfunctions. First, wc construct a for- 
mal eigenfunction of the operator T'^{w, Dy,) admitting an asymptotic expansion 
in the form of a formal asymptotic series in powers of h^^^ 

00 

1=0 

with the corresponding formal eigenvalue 

00 

such that 

T^{w,D^)u'' -X'^u'' = Q 
in the sense of asymptotic series in powers of h}/'^. 

The first terms. Looking at the coefficient of hP, we obtain: 

To(x2,i5.J«W=AW«W. 

Thus, we have 

(2.10) AW=Af = (2fc + l)6o, u^°\x2,^) = ^Mx2)xa{0, k&n, 
where xo is some function, which will be determined later, and c is some constant. 
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Looking at the coefficient of /i^/^, we obtain: 

The orthogonaUty condition imphes that 

A(i) = 0. 

Under this condition, we get 

(2.11) ««(^2,0 = ^V'fe(a;2)xi(0, 

where Xi is some function, which will be determined later. 
Next, the cancelation of the coefficient of gives: 

(2.12) {To{x2, D,,) - aW)«(2) = A(2)wW - T^iw, D^)u^°l 
The orthogonality condition implies that 

A^'^X0(O - TT^ I T2{w,D^)u^''^i,k{X2)dX2 = 0. 

W^kW J 

Lemma 2.5. For any function u of the form u{x2,£,) = p]^'*/'/c(a;2)x(?)j '''e have 

||^-|| j T2{'w,Dyj)u{w)ipk{x2)dx2 ='Hfcx(0. 
where Hk is the harmonic oscillator: 

Hk =(2fc + + {2k + l)ai-^e^ 

+ ^(2fc' + 2A: + 1) (^ai + h + h}oR{xo)^ - -^R{xo). 

Proof. We have 

Dx2'^k = -i^y^ (2kipk-i - bl^^X2ipkj , 
Dl^jpk = bo{2k + l-boxl)jl>k, 

and 

Dl^tPk = 2iblx2i^k + bo{2k + 1 - box^D^^TPk- 

We also have 

2xHk =Hk+i+2kHk-i, 

Ax'^Hk =Hk+2 + (4fc + 2)Hk + 4k{k ~ l)Hk-2, 

Sx^Hk =Hk+3 + (6fc + 6)fffc+i + 12k^Hk-i + 8k{k ~ l)(fc - 2)Hk-3 
16x'^Hk{x) =Hk+4{x) + {8k + l2)Hk+2{x) + 12(2fc2 + 2fc + l)Hk{x) 

+ 16(2A:2 - 3A: + l)kHk-2{x) + l&k{k - l){k - 2){k - i)Hk-i{x), 
that implies that 

1 1 

;{x2i}k-i,i>k) 



\Hk\r-^''"'-''^'"-2by^' 

^ :{xltpk-i,tpk) =-^k, 



j^ixli^k,^} =2^(2fc + l), 
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^ {xt^k,i^k) =^{2k'' + 2k + l). 



Next, we have 



= 2^' 

^ 'n2 „/, „/, \ _ ''0 //oi. II 



:{Di^i>k,i>k) =p^((2fc + l-6oar^)Vfe,Vfe) 
=^(2fc + l)6o, 



{x2D^^i)k,i)k) = iiu 112 (((^^ + 1)^2 - ''oa;4)V'fe, V'fe) 



ll^felP^ ' ' ll^fcll 

=^(2fc2 +2A:- 1), 

First, remark that 



=-(2fc2+2fc + l)62. 



^ ]|^^ / (~^*'^^^°^^'^'^^ ~ ^iR(xo)bQX2D^ u{w)'il)k{x2)dx2 = 0, 

since the operator —^iR{xo)-^Dx2 — ■^iR{xo)boX2D(^ is odd in the X2 variable. 
Next, we have 

^lA =^^{xli,k,i>k)mD^)x{0 
={2k + l) (^^^e+PDl^ X(0- 



Now we consider 

^2A=T7^-n- / bo{X2M{w,Dyj) + M{w,Dyj)x2)u{w)ll)k{X2)dX2. 
W^kW J 

The operators Mi{w, D^) and M2{w, D^) are odd in the X2 variable, so we obtain 

•^2^=77^77 / bo{X2Mi{w,D^) + Mx{w,Dw)x2)u{w)'ll)k{X2)dX2 

\\Uk\\ J 

bo {X2M2{W, Dyj) + M2(W, D^^)X2) u{w)lpk{X2)dX2 



Hk 



iC^xi + ToADI^ + ToX2Dlx2 ) Vfc, Vfc)x(0 



\\Hkr\i-^^^ br^^^^^ bi 

~ J^{(^^X2D.,,+2^D.,,X2 ) i^k.i^k)iD^x{0 
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26o Y 2/3 ^ , ./3' 



Thus we arrive at 
Next, we consider 

By (|2.5p . we have 

9ih + -i?e + ^i?..) = IrM{-d, + ^D^^f. 

Therefore, 

It suffices to consider the terms, which are even with respect to X2'- 
S3X - ^^^(^0) J {Dl + ^^Dl;)Dl^u{w)Mx2)dx2. 



Thus, we obtain that 

I 

6||Ffc 

I 



(ibl\\Hk\ 

(2 
I 



:l-{2k + l)boRixo)Dlx{0 



-(2fc2 + 2fc+I)i?(a;o)x(e). 
o 

Next, we consider 

S4X = ^"1"^^ J ^^^^ (^^ ^ V' ^^^ ^ V^'^'^ (^Dx^{boX2)+boX2DrcJju{w)ipk{x2)dx2. 
By (|2.5p . we have 

'^4A = ^j^^-jyi?(a;o) y" ^2^2 + (^^* ^ ^"^^^) (2^o^2£'2;2 " ibo) u{w)'iljk{x2)dx2. 
It suffices to consider the terms, which are even with respect to X2: 

'54A =^j^^--j|i?(xo) J {x2Dx.^ - S,D{}{2x2Dx^ - i^u{w)tpkix2)dx2 



3||7?fc 



2 



iR{xo){x2D^^'ipk,'ipk)x{0 



2\\Hk\ 

-^l-^R{xo){{2x2D,,-t)i^k,^k)^D^x{0 
^(2fc2 + 2fc + l)i?(xo) + ^R{xo)] xiO- 
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By (^31) . we have 

= rr II -^(^o) / (2^2 + ■^ybQX2u{w)^/jk{x2)dx2. 

It suffices to consider the terms, which are even with respect to X2'- 



= cll u_ II R{xa) I {x\ + -T0^\x\'^{w)'^k{x2)dx2 



Next, consider 



By (P^ . we have 



y"r(i) + ^, -I?^ + ^Dj,^ Dr,^u{w)tpkix2)dx2- 



^ 3||ff^,|| ^('^o) y ^2^2 + Drc^u{w)ipk{x2)dx2. 
It suffices to consider the terms, which are even with respect to X2'- 
'^eA = Qiirj ||2 ^('^o)(a:2-Px2V'fc>fc)x(C) 

= ^Ir{xo)x{0- 



Finally, we take 

i 



5j\ = 



By (P^ . we have 



y"r(;^) + y,-D^ + ^Dx2^ hoX2u{w)ll)k{x2)dX2. 



SiX = 3||ff^|| ^(^o) y" (^"^^ ^ V^^^^ hQX2u{w)'4>k{x2)dx2. 
It suffices to consider the terms, which are even with respect to X2'. 

= -^i?(2;o)x(0- 



We conclude 

1 



T2{w, Dw)u{w)tpkix2)dx2 



\\Hk\ 

=doX + SiX + S2X + S3X + S4X + S^X + SqX + SrX 
= ^(2fc' + 2k+l)(a + p+ hoRixo)) xiO 
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+ {2k + l)[l3 + ^boR{xo) ) DixiO 



12 

= ^(2fc2 + 2fc + 1) (ai + /3i + ^boR{xo)^ xiO 
+ {2k + l)PiDlx{0 + (2fc + Ihi^exiO - \r{xo)x{0- 



□ 



Thus, we obtain that 



X^'^ = A',} = '^.k, J, fee 



where Vjk is an eigenvalue of the harmonic oscillator "Hfei 
iy,k - (2j + l)(2fc + l)(ai/3i)i/2l 

+ 7^(2A:2 + 2fc + 1) (ai + /3i) + ^(fc^ + k)R{xo), j, fc G N, 

and 

where ^jfc is the normalized eigenfunction of Hk associated to the eigenvalue Vjk- 
Moreover, we conclude that u^^^ is a solution of (|2.12p . which can be written as 

where 0^^-' is a solution of ()2.12p . satisfying the condition 

C^^^Hx2,OMx2)dX2^0, 

and X2 will be determined later. 

Now the cancelation of the coefficient of /i'^/^ gives: 

(2.13) (To(a;2,i?x.)-A("')u(3) 

The orthogonality condition for (|2.13p is written as 

(2.14) A(3)xo - J niw, D^)u^''^M^2) dx2 + A(2)xi - HkXi = 0. 

Under this assumption, we obtain that u*^^^ is a solution of (|2.13p . which can be 
written as 

7.(3) = 0(3)(X2,O + ^^fc(x2)X3(e), 
ll-Wfcll 

where (j)^^^ is a solution of (|2.13p . satisfying the condition 

0(3)(x2,OV'fc(a;2)rfa;2 =0, 
and X3 will be determined later. 
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The equation p.l4[) has a solution if and only if 

A(3) = JW^JJ T3{w,D^)u^°\w)Mx2)^jk{0dx2d^, 

which allow us to find A^'^^. Under this condition, there exists a unique solution xi 
of (|2.14p , orthogonal to xo 

The iteration procedure. Suppose that the coefficients of h^^^ equal zero for 
i — 0, . . . ,n — 1, n > 3. Then we know the coefficients A'-^'' for £ = 0, . . . , n — 1. We 
also know that u'^' for ^ = 0, . . . , n — 1 can be written as 

\\J^k\\ 

where (f)'-^\i = 0, — 1, are some known functions in iS(M'^), satisfying the 

condition 

and xe € §(R) are known for £ = 0, . . . , n — 3, x« -L Xo- 
The cancelation of the coefficient of ft."/^ gives: 

(2.15) (To(a;2,i?..)-A("')u(") =A(")7.W-r„(7«, 

n-l 

The orthogonality condition for (|2.15p is written as 

(2.16) a(")xo- J Tn{w,D^)u'^°'>i;k(x2)dx2 

n-l „ 

+ ^(AWx«-£-^F^ / Teiw,D,,)u^''-'^k{x2)dx2) 



\Hk\ 



1 



Te{w, D^)cj)(^"'''>M^2) dX2 + \^^hn-2 - HfeXn- 



Under this assumption, we obtain that u^"^ is a solution of (12.151) . which can be 
written as 

= 0(")(^2,O + Tr^Mx2)Xn{0, 
W^kW 

where (j)^"'^ is a solution of (I2.15p . satisfying the condition 

^^"\x2,OMx2)dX2^0, 

and Xn will be determined later. 

The orthogonality condition for (|2.16l) allows us to find A*^"^ . Under this condi- 
tion, there exists a unique solution Xn-2 of (|2.16p . orthogonal to xo- 

Thus, for any j G N and k E N, we have constructed an approximate eigenfunc- 
tion u^j. of the operator T^{w, D^) admitting an asymptotic expansion in the form 
of a formal asymptotic series in powers of /i^/^ 

oo 
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such that 

(2.17) ufjix,,0 = Tr^M^2)^M0- 

ll-Wfell 

with the corresponding approximate eigenvalue 

oo 

For any N G N, consider 

N N 

Then we have 

The constructed functions w^j.^^^ have sufficient decay properties. Therefore, by 
changing back to the original coordinates and muhiplying by a fixed cut-off function, 
we obtain the desired functions (/'^j.^, which satisfy (|2.2p with ^j,k.i — ^f^- 

The system {uf^} is an orthonormal system. Since each change of variables, 
which we use, is unitary, this implies the condition (|2.ip . 

3. Lower bounds 

In this section, we will prove the lower bound in Theorem ll.il First, we recall a 
general lower bound due to Montgomery [37]. Suppose that J7 is a domain in M. 
Then, for any u €E C^{U), the following estimate holds: 



(3.1) \\{ihd + A)u\\l> 



b\u\'^dxg 



This fact is an immediate consequence of a Weitzenbock-Bochner type identity. 

From (|3.ip . it follows that we can restrict our considerations by any sufficiently 
small neighborhood ft oi xq. Denote by the Dirichlet realization of the operator 
in L^{n,dxg). 

The estimate p.ip implies that 

ThH'^ + {l-Th)hb<H'l,, 0<T<h-\ 
Taking r — h~^/^, we obtain 

fll/2(^jfh _ j^^ ^ ^1/2^) < ^ < /i< 1. 

Consider the Dirichlet realization of the operator i/'* — hb + h^/^{b — bo) in 
L^(f7, dxg). Then we have 

(3.2) hbo + h^^^Xj{P^) <Xj{H^). 

Therefore, the desired lower bound for Xj{H^) is an immediate consequence of the 
following theorem. 

Theorem 3.1. For any j e N, there exist Cj > and hj > such that 

,2 n 



2d^. 
bo ^^2bo 



C-jh^^"'^, he{0,hj 
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To prove Theorem lS-H we will follow the lines of the proof of 11, Theorem 7.4]. 
First, we observe that the upper bound in Theorem 11.11 and p.2p imply an upper 
bound for XjiPji,): 

-J + 



(3.3) Xj{P^) < /i3/2 



bo 2bo\ 

For any eigenvalue Xj{P!^), denote by u'}^^ an associated eigenfunction. By a 
straightforward repetition of the arguments of [H] , we can easily show the following 
analogue of Lemmas 7.10 and 7.11 in llj^. 

Lemma 3.2. For any j € N and any real k > 0, we have 
For any j e N, any a = 1, 2 and any k > 0, we have 

Take normal local coordinates near xo such that xo corresponds to (0,0) and, in 
a neighborhood of xq, 

b{X) = bo + aix^ + Piy^ + 0(|X|3). 

So we have 

5ii(X) = l + 0(|Xn, g,,{X) = 0{\X\'), g22{X) = l + 0{\X\'). 
We can take a magnetic potential A such that 

Ai{X) = -^boy + 0{\Xf}, A^iX) = ^box + 0{\Xf}. 
Let us introduce 

b2iX) = -X - Hess 6(0) X. 

Thus, we have 

b{X) = bo + b2{X) + 0{\X\^). 

We have 

Pd^ 9""iX)^'L^'p + ^h ^ V^Vl-hb{X) + h'/\b{X)~bo), 

l<a,f}<2 l<a<2 

so its quadratic form is given by 

-h f b{X)\u{X)\^.Mx)dX + h^/^ f {b{X) - bo)\u{X)\''^^dX. 
Jn Jn 

Note that 



Pd > 0- 

Now we move the operator into the Hilbert space L^{il,dX), using the 
unitary change of variables v — \g(X')\^l'^u. For the corresponding operator 

Pi = \g{X)\^I^Pl\g{X)\-^l^ 



^There are a few inaccuracies in concerning Lemma 7.11. For the erratum, see 

Ihttp:/ /www. math. u-psud.fr/~hclffer/erratuml64II.pdf 
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in i^(ri, dX), wc obtain 



{P!^v,v)=[ 9"^{X)(v'^„-lih\giX)\-'-^\giX)\)v{X)x 

•^^ l<a,l3<2 ^ a J 

X (^\--^h\g{X)\-^^\g{X)^v{X)dX 

-hi h{X)\v{X)\^dX^-h}l'^ [ {b{X)-bo)\v{X)fdX. 
Jn Jn 



Put 



q{v) 



E 

l<a<2 



V'^^-iih\g{X)\-'-^\giX)\]v{X) 



dX. 



Then we have 
(3.4) \{P^v,v)-q{v 

< 



■1^ 1^2 ^ ^ 



h\g{X)\-^^\g{X)\ ) ^(X) 



Ci/i / \v{X)\''dX + C2h^/^ [ \X\^\v{X)\^dX. 



Consider the Dirichlct rcahzation P^-'i^^, jj of the operator 



d__l 



boy 



in the space L^(Q, rfX). So its quadratic form is given by 

{P'}i,^v,v) = qf"^\v)-hho f \v{X)\'dX + h'/' f b2{X)\v{X)\'dX, 

Jn Jn 



d 1 



dX. 



where 

So we have 

(3.5) \{P^v,v) - {P}i,,v,v)\ < \q{v) - qf''^\v)\ 

+ V / \X\^\\/lv{X)\^dX + h [ \X\^\v{X)\^dX + h'/^ [ \X\^\v{X)\^dX. 
^ Jn Jn Jn 

Finally, we have 

(3.6) \q{v) - qf^'^'iv)] 



<C{q{v)f" 



\x\Mx)?dx 



1/2 



\x\'Hx)\^dx 



1/2- 



For a fixed j G N, consider the subspacc V^'^ of L'^{il,dxg), generated by all 
eigenfunctions of the operator Pj^ associated to the eigenvalue Af(P^) with i = 
0, 1, . . Thus, V^'^ is a (j + l)-diniensional space such that 

(3.7) {Phuh,UH) < XjiP^))\\uH\\h^n,a,^y uneV^'^. 
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Moreover, by Lemma for any real fc > 0, there exists Cfe > such that, for any 

(3.8) \\\X\'uhh^n,d.,) < Cfc(/i*^/2 + /i(3'=+i)/«)h„|U.(n,d.,), 
for any a — 1,2 and any fc > 0, 

(3.9) WlXl'^^'lu^h^nM.,) < Ck{h'^''+'^/' + h(''^+'y^)\\u,,h^n,d.,). 

By p.Sp and p.9p . for any real fc > 0, there exists Ck > such that, for any 
Vh £ L^{n,dX) of the form Vh = \g{X)\^/W with Uh € V^'^ , 

(3.10) \\\X\'^Vk]\L^in,dX) < Ckih""/^ + h^''^+^y^)\\vhU^nMX), 
for any a — 1,2 and for any fc > 0. 

\ 1/2 



(3.11) 



1^1 



2k 



d 1 



1^1 



2fe 



d 1 



2 \ 1/2 



<C,(/.(^-+l)/2 + /i(3'^-+5)/8)||„^||^,^^^^^^. 

The estimates p.lOp and p. lip allow us show for any Vh G L?'{Vl,dX) of the 
form Vh = \g{X)\^/^Uh with uh e F'^^, first, using dSJ]), dSS]) and dSH), that, 

q{vh) < C/i3/2|| 

^'i|lL2(a,dx)' 

next, using p.6p . that 

and finally, using p. 51) . that 

(3.12) (P^t"/.,"/.) < (A,(P5) + Q/ii^/«)|| 

Let X be a function from C^(M") such that suppx C fl and x = 1 in a neigh- 
borhood of zero. By p.lOp and (|3.11l) . it follows that, for any fc G N there exists 
Cfc > such that, for any Vh G L^{fl,dX) of the form Vh = \g{X)\^/'^Uh with 

(3.13) 11(1 - x)vh\\LHn,dx) + \\-Q^Vh\\L^{n,dx) + W^v^W L'^in.dx) 



< 



C'kh''\\vh\\L^(Q.dX)- 



and 



(3.14) 



a-xix)) 

+ 1 f (i-xix)) 



d 1 



1/2 



dX 



ih^^ + ^box ) Vh{X) 



2 \ 1/2 
dX 



< Ckh''\\vh\\L^{n,dx)- 
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Using (|3.13|) and p.l4p . it is easy to check that, for any A: > 0, there exists 
Cfe > such that 



Consider the self-adjoint reahzation of the operator Pj/^j in L"^ {M? , dX) . We 
will the same notation P^i^^ for this operator. Consider the (j + l)-dimensional 
subspace W'''^ of C^iR"^), which consists of aU functions Wh e C^{R^) of the 
form Wh = x\9{X)\^^W with u,, € V''''^ . By (jXT^ and (jXTSj) . it follows that, for 

any Wh G W'^'^ , 

{P^l,,WH,Wu)<{HPh) + C,h^''l^)\\ 

^'i|li2(o,<ix)- 

By the mini-max principle, this immediately implies that, for any j > 0, there 
exists Cj > such that, for j-th eigenvalue ^j{Pfiat) of ^fiaf have 

(3.16) A,(P;,,J<A,(Fi5)+Q/ii5/«. 

It remains to recall that the eigenvalues of the Schrodinger operator with constant 
magnetic field and positive quadratic potential in M" can be computed explicitly. 
More precisely (see, for instance, [26l Theorem 2.2]), the eigenvalues of the operator 



d 1, 



d 1, 



dy 2 



are given by 
where 



{2ni + l)si + (2^2 + 1)S2, ni, 712 e N, 



S2 



1 
1 

71 



1/2 



1/2 



and 



tx^T^rK, dK = detK. 

flaV 

A„i„2 = (S'T-l + l).5l + (2^2 + 1)S2 - ^^0, ^1, 712 G N, 



Applying this formula to the operator P^i^^^ , we obtain that its eigenvalues have the 
form: 

where 



si 



h^'h + bl- [{h^/h + blf - Ah^'^d]"^ 
V2 L 

d^'H^^h'/^ + Oih^), 



1/2 



and 



32 



V2 



h^'h + bl + [{h^'h + blf -Ah^'M^/' 



.1/2 



hbo + ^tbo^h^^^ + 0{h^). 



Thus, we obtain 

Anin2 — 27i2/l6o 



■ 2dl/2 t ^2 

oq oo 26o 



/l'^/2 + 0(/l2)^ 711,772 eN. 
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For jth eigenvalue ^j{Pfiat) of ^fiaf '^^ obtain 

(3.17) XjiP^lat) = 

Conibining p.l6p and p.l7p . we immediately complete the proof of Theorem 13. II 

Proof of Theorem [TE Fix j e N. By Theorem [TTTl there exist C > and /iq > 
such that, for any h £ (0, ho], 

I, fl Spec{H'^) = {X,{H^)}, 

where 

Ij = (^hbo + h^ 

On the other hand, by Corollary I2.2[ for any natural iV, there exist C" > and 
/ig > such that, for any h G (0, /iq], 

dist(^^o^,Spec(i7'')) < C'h^. 
Without loss of generality, we can assume that, for any h e (0, min(/io, /iq)], 

{^Ji%J, - c'h'^,fi%j^ + c'h"^) n /, - 0, V£ ^ J. 

Hence, for any h e (0, min(/io, /iq)], Xj{H'^) is the point of Spec(i?'*), closest to 
t^%N- It follows that 

\X,{H'') - < C'h^, h&{0,min{ho,h'o)], 

that proves (jl.ip with ajj = fJ-j,o,i- D 

4. Periodic case and spectral gaps 

In this Section, we apply the results of Section [5] to the problem of existence of 
gaps in the spectrum of a periodic magnetic Schrodinger operator. Some related 
results on spectral gaps for periodic magnetic Schrodinger operators can be found 
in[ll[Sl[IillIll[I71[IHl[ini[2ni|2Il[ll|ll[ig (see also the references therein) . 

Let M be a two-dimensional noncompact oriented manifold of dimension n > 2 
equipped with a properly discontinuous action of a finitely generated, discrete group 
r such that M/T is compact. Suppose that iJ^(M, M) = 0, i.e. any closed 1-form 
on M is exact. Let g be a F-invariant Ricmannian metric and B a real-valued 
F-invariant closed 2-form on M. Assume that B is exact and choose a real-valued 
1-form A on M such that dA = B. Write B = bdxg, where b G C°°{M) and dxg is 
the Riemannian volume form. Let 

&o — min b(x). 

xeM 

Assume that there exist a (connected) fundamental domain T and a constant eo > 
such that 

b{x) >bQ + eo, X <E dT . 

We will consider the magnetic Schrodinger operator as an unbounded self- 
adjoint operator in the Hilbert space L'^{M). Using the results of [5], one can 
immediately derive from Theorem 12.11 the following result on existence of gaps in 
the spectrum of in the semiclassical limit. 



2d^ . 

~br^^2b'o 



h^/^ + 0{h^). 



2d^ . 



Ch'^/\hbo + h^ 



2 



2d^ . a_ 



Ch^'^ 



22 



BERNARD HELFFER AND YURI A. KORDYUKOV 



We will use the above notation 

t = Tr Qhbss 5(2:0)^ , d = dot QHess5(a;o)^ . 
For any fc S N, put 

+ 1 

Ck = {2k + 1)— + (2fc2 + 2k + 1)— + -(fc2 + k)R{x„). 
oo 2bo I 

Theorem 4.1. Assume that &o > and there exist xq ^ J- and C > such that 
for all X in some neighborhood of xq the estimates hold: 

C-^ d{x, xo)^ < b{x) -bo<C d{x, xo)^ . 

Then, for any natural k and N , there exist Ck,N > Cfc 0,'rid hk,N > such that the 
spectrum of in the interval 

[{2k + l)hbo + h^Ck, {2k + l)hbo + h^Ck.N] 

has at least N gaps for any h G (0, ft.fe,Ar]. 
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